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We introduce tangent space methods for projected entangled-pair states (PEPS) that provide 
direct access to the low-energy sector of strongly-correlated two-dimensional quantum systems. 
More specifically, we construct a variational ansatz for elementary excitations on top of PEPS 
ground states that allows for computing gaps, dispersion relations and spectral weights directly in 
the thermodynamic limit. Solving the corresponding variational problem requires the evaluation 
of momentum transformed two-point and three-point correlation functions on a PEPS background, 
which we can compute efficiently by introducing a new contraction scheme. As an application 
we study the spectral properties of the magnons of the Affieck-Kennedy-Lieb-Tasaki model on the 
square lattice and the anyonic excitatons in a perturbed version of Kitaev’s toric code. 


In the last decades, low-dimensional quantum systems 
have been at the forefront of both theoretical and ex¬ 
perimental physics. With the reduced dimensionality al¬ 
lowing for stronger quantum correlations, these systems 
show an extreme variety of exotic phenomena but are 
notoriously difficult to simulate. In particular, the low- 
energy physics of such systems typically cannot be ap¬ 
prehended by perturbing around some free or otherwise 
exactly solvable point. 

For one-dimensional quantum spin chains, the advent 
of the density matrix renormalization group (DMRG)P 
has proven revolutionary. DMRG has developed into 
the de facto standard method for reliably and efficiently 
probing the low-energy behaviour of strongly correlated 
systems in one dimension. A better understanding of 
this success was realized through the identification of the 
variational class over wh ich DMRG optimizes as matrix 
product states (MPS;P®. More specifically, it was un¬ 
derstood how the entanglement structure of MPS allows 
for a natural parametrization of the low-energy states of 
gapped, local one-dimensional quantum systemj^. 

The characterizing feature of the entanglement struc¬ 
ture in gapped local quantum swtems is the observed 
area law for entanglement entropy^. With this insight the 
natural extension of MPS to higher dimensions is gi ven 
by the class of projected entangled-pair states (PEPS)^. 
This set of states has shown to capture the low-energy 
physics of several interesting systems in two dimensions 
and is competing with more establishe d m ethods in de¬ 
termining their ground-state propertiej^ES! Yet, in con¬ 
trast to the one-dimensional case, the PEPS simulation 
of two-dimensional systems is computationally challeng¬ 
ing, thus making the development of new algorithms 
highly desirable. Eurthermore, existing PEPS algorithms 
focus exclusively on capturing ground state wave func¬ 
tions, whereas experimentally relevant low-energy prop¬ 
erties such as excitation spectra remain inaccessible. 

In this paper we initiate a new set of methods based on 
the tangent space of the PEPS manifold of states, a con¬ 
cept which has proven extremely versatile in the context 
of matrix product stateJ^^. We apply this methodology 
for constructing a variational ansatz for elementary ex¬ 


citations on top of PEPS ground states. Combined with 
the versatility of the PEPS as ground state ansatz, our 
method provides a flexible and mostly unbiased approach 
(aside from entanglement considerations in the PEPS 
ansatz) towards extracting the low-energy spectrum of 
a strongly correlated two-dimensional quantum system. 
We explain the main ingredients for successfully apply¬ 
ing the variational principle, for which we need to intro¬ 
duce a new PEPS contraction scheme. As an application 
we study the excitation spectrum of the two-dimensional 
Affieck-Kennedy-Lieb-Tasaki (AKLT) model on a square 
lattice, and of a perturbed version of the toric code 
model, where we show the ability to access single topo¬ 
logical excitations (anyons) directly. 

Elementary excitations. We start from one of the cen¬ 
tral insights of condensed-matter physics that the low- 
energy properties of quantum systems can be understood 
in terms of elementary excitations or quasi-particleJ^^. 
Experimental observables such as dynamical correlation 
functions, real-time evolution or low-temperature prop¬ 
erties can be understood from this picture of weakly- 
interacting particles on a non-trivial vacuum state. This 
insight can be materialized within the PEPS formalism 
by first identifying the PEPS class of states as a vari¬ 
ational manifold, embedded in the full Hilbert space, 
that captures the ground state structure of local Hamil¬ 
tonians. Analoguous to e.g. the manifold of Slater 
determinants for the Hartree-Fock approximation of a 
fermionic ground stat^i^, or the manifold of matrix prod¬ 
uct states^ for ground states of quantum spin chains, 
the excitations on top of such a ground state are ob¬ 
tained as linear perturbations living in the tangent space 
to the manifold rather than the manifold itself. This gives 
rise to a set of met hods k nown as post- Hartree-Fock^ or 
post-MPS method^iiE6|^ respectively. Here, we develop 
“post-PEPS” methods by constructing elementary exci¬ 
tations within the tangent space of PEPS ground states. 

PEPS and variational excitations. Gonsider a two- 
dimensional quantum lattice system in the thermody¬ 
namic limit. A translation-invariant PEP^ can be 
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parametrized by a single five-legged tensor r di 

d 

|VI/(A)) = 5] C2({^«}) |{s}) 

14=1 

where C 2 (...) denotes the contraction of an infinite two- 
dimensional network of A tensors, which is more conve¬ 
niently represented pictorially as 


MA)) = 



( 1 ) 


Here, we have represented the tensor A as 

u 

■^u.l.d.r , /T 

d ,s 

The dimension of the four virtual legs {u, r,l,d) is called 
the bond dimension D of the PEPS and serves as a re¬ 
finement parameter of the variational class, whereas the 
leg s has the local physical dimension. In diagrams such 
as in Eq. Q whenever two legs are connected the two 
corresponding indices are contracted (i.e. identified and 
summed over), such that |4'(H)) is obtained by contract¬ 
ing virtual legs of neighbouring tensorJi^. 

Computing the norm or a local expectation value of a 
PEPS involves the contraction of an infinite tensor net¬ 
work. Indeed, denoting the “double layer” tensor a as 

O-uu' ,rr' ,dd' ,11' — ^ ® = 

s 

we can represent the norm of the state in “top view” as 





This infin ite contraction cannot be performed 
exactlj®^ but, by relying on MPS-style algorithms, a 
precise approximation can be computed efficiently, as 
discussed below. With a slight modification, the same 
algorithms can be used to compute the expectation value 
of local observables, such as the energy density, and 
therefore suffice to implement the variational principle. 

Suppose we have found a translation-invariant PEPS 
representation for the ground state of a local Hamiltonian 
H. The tangent space ansatz for elementary excitations 
is given by 


m,n 



( 3 ) 


where the red tensors correspond to those of the optimal 
ground state and only the blue circle represents a new 


tensor B at site (m,n). This ansatz thus represents the 
momentum superposition of a local perturbation, which 
has a finite width (determined by the PEPS bond di¬ 
mension) by acting via the virtual de grees of freedom. It 
encompasses the Feynman-Bijl ansatd^filSU and its valid¬ 
ity for gapped excitations can be rigorously motivatecP^. 

Once we have defined the variational subspace, we can 
minimize the energy in order to find the best approxima¬ 
tion to the true excitation. Since the subspace is linear, 
this gives rise to the Rayleigh-Ritz problem 

(4) 

where and are the effective Hamiltonian 

and norm matrices, which can be interpreted as carry¬ 
ing the one-particle dynamics on top of a PEPS back¬ 
ground. The lowest eigenvalue ui corresponds to the low¬ 
est excitation energy at momentum {Kx,Ky). By repeat¬ 
ing this procedure for different momenta, we can obtain 
the dispersion relation for all one-particle excitations in 
the system. Moreover, with an expression for the wave 
functions, we can straightforwardly compute their spec¬ 
tral weights. 

Effective environments. The matrix elements of the 
effective Hamiltonian and normalization matrices contain 
Fourier transformed two- and three-point functions and 
cannot be computed efficiently using the current PEPS 
contraction schemes for local expectation values. PEPS 
contraction schemes try to capture the effect of all the 
surrounding tensors in an effective environment using an 
approximate model. Depending on the situation, the 
most efficient schemes eith er model the environment as 
a matrix product stat^^EH or make use of the idea of a 
corner transfer matrixP^t^. For the problem at hand, 
we combine ideas from both schemes to construct a new 
contraction method. 

Let us first review the most straightforward algorithm 
for contracting the network Eq. ([^, which consists of 
running through the network in one direction, sequen¬ 
tially applying a row of tensors, and approximating the 
boundary as an MPS in every step. In the infinite case, 
where we have an infinite number of rows, this algorithm 
boils down to finding a translation invariant MPS repre¬ 
sentation of the fixed point of the linear transfer matrix 
T, which can be pictorially represented as 



Once we have found this fixed point and rescaled the 
PEPS tensor such that the largest eigenvalue of T is 
1 (corresponding to a normalized PEPS), the computa¬ 
tion of local expectation values is analoguous to well- 
established and efficient MPS algorithms. One can also 
construct approximate eigenstates for the next eigenval¬ 
ues of the transfer matrix T and infer from it the qual ita- 
tive features of the physical excitation spectruni2SEZl_ 
particular, the second largest eigenvalue Xf' is related to 
the correlation length of the PEPS as ^ = — 1/log |A^^|. 
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The computation of general n-point functions are more 
difficult (unless the n points are aligned along one of the 
main axes), but are a necessary ingredient in the con¬ 
struction of the effective matrices in Eq. We there¬ 
fore introduce a new contraction scheme: the central idea 
is that, instead of contracting Eq. ([^ top-down linearly, 
we run through the lattice diagonally by sequentially ap¬ 
plying a corner-shaped transfer matrix C 


C = 



We can again approximate the “fixed point” of this op¬ 
erator as a boundary MPS. One can verify that the only 
difference with the linear case is the presence of the cor¬ 
ner, which can be captured by the insertion of a corner 
matrix 5, which acts on the virtual level of the uniform 
MPS and is depicted with a diamond shape below. The 
fixed point equation then reduces to a linear equation in 
S so that the optimal solution is easily founcP^. 

Using the fixed points of all four corner transfer matri¬ 
ces, we can construct an effective environment consisting 
of different one-dimensional ” channels”, for which we can 
compute expectation values exactly using standard MPS 
technique's. More specihcally, with the channel struc¬ 
ture we can compute two- and three-point functions for 
operators at arbitrary relative positions by contracting 
networks such as 



where the colored tensors indicate where the operators 
are located. In the Supplemental Material we show how 
to compute fourier transforms of these types of diagrams 
directly, thus enabling the evaluation of e.g. static struc¬ 
ture factors and the matrix elements featuring in Eq. Q . 

Magnons in the AKLT model. As a first application, 
we study the two-dimensional AKLT modeP^ on the 
square lattice. The one-dimensional AKLT model was 
first introduced to establish the existence of a gap in 
a rotationally invariant spin-1 chain, in support of the 
Haldane conjectur^^. The construction can be straight¬ 
forwardly extended to two-dimensional systems, but the 



FIG. 1. Single mode approximation for the dispersion relation 
of the two-dimensional AKLT model on a square lattice. The 
minimum of the triplet dispersion relation is at momentnm 

^x,y ~ • 


existence of a gap - although highly expected - can no 
longer be established rigorously. On the square lattice 
the AKLT state can be represented as a PEPS with bond 
dimension = 2, and is the unique ground state of 

^^AKLT = ^ E + ^h% 

<b) 

with hij = Si ■ Sj the spin-2 Heisenberg interaction. 

The elementary excitations in this model are expected 
to be (triply-degenerate) magnons. In a first step, we can 
target the magnons with the single-mode approximation 
(SMA), which is known to reproduce the one-particle 
dispersion relation qualitatively in the one-dimensional 
casd^. In the present context the SMA wave function is 
given by (a = x, y, z) 

m,n 

The norm of these states ) is equal to the 

static structure factor and can be computed with our new 
contraction schem^^ while the energy expectation value 
reduces to an easy PEPS contraction. The dispersion 
relation is shown in Fig. the spectrum consists of an 
elementary triplet with its minimum at Kx,y = Tp^and a 
gap Asma = 0.0199. 

Next we can determine the excitations variationally 
with the excitation ansatz [Eq. ([^]. To improve on the 
SMA result, we enlarge our variational subspace by per¬ 
turbing the ground state on a larger region - a procedure 
that is bound to converge exponentially fast to the cor¬ 
rect wave functiorP^. By introducing a B tensor on a 
block of two by two sites, we are able to estimate the 
gap at Avar = O.OMiPI, in excellent agreement with the 
value 2A = 0.03 in Ref. [MJ obtained through a com¬ 
putation of the magnetization curve. In addition, we 
can compute the characteristic velocity through the sec¬ 
ond derivative of the dispersion relation in its minimum 
and we obtain Uvar = 0.04115. Finally, in Table we 
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p 

spectral weight w 

% of sum rule 

O.OOOOtt 

31.222 

99.57 

0.06 827r 

26.058 

99.35 

0.13647r 

17.198 

98.84 

0.20 457r 

10.672 

96.59 

0.22 737r 

8.9968 

94.26 

0.25 007r 

7.3452 

88.75 


TABLE I. Spectral weight w = K'^“l l^(^))l^ ^he 

elementary magnon triplet as a function of momentum = 
Ky = TV + pIV2, and the percentage of the integrated spectral 
function that is saturated in the one-particle sector (we can 
compute the integrated spectral function exactly through the 
sum rulj^ and the static structure factor, see Supplemental 
Material). We observe that the magnon contains nearly all 
spectral weight in the minimum, and that this percentage 
drops when going away from this point. 

have listed the spectral weight of the elementary magnon 
around its minimum. 

Anyons in the perturbed toric code. As a second ex¬ 
ample we study the anyonic excitations in the toric 
cod^^, the easiest example of a model exhibiting topo¬ 
logical order and anyonic excitations (fluxes and charges). 
At the toric code fixed point, the anyons have a flat dis¬ 
persion relation since all Hamiltonian terms commute. 
To generate non-trivial dynamics, we perturb the toric 
code state corresponding to a non-unitary operator 

fi = exp Q(/3rr< + /dzcr-)^ . 
acting on every site i. The associated Hamiltonian is 
^FTc +y^Hp 

S p 

where the filtered star operators are given by 



and similar for the plaquette operator Hp with for 
i £ p. In first order, the filtering operation is equivalent 
to applying a magnetic field to the toric code. In fact, 
it has been shown in Refs. [27] and 1371 that the phase 
diagram is qualitatively similar. 

In a number of recent works it has been established 
that PEPS can provide a natural description of topo¬ 
logical phases, where the topological order of the global 
state is reflected by a sy mmet ry of the local PEPS ten¬ 
sor A on the virtual leveP^^^. A variational ansatz for 
the complete set of elementary excitations in the different 
anyon sectors is obtained by attaching a half-infinite vir¬ 
tual matrix product operator string to the local B tensor 
in the excitation ansatz [Eq. ©]• The topological sector 
is encoded in the type of string and the virtual symme¬ 
try representation of the local tensor. In the case of the 



FIG. 2. The dispersion relation of the toric code in the flux 
(blue-green) and charge (red-yellow) sector for filtering pa¬ 
rameters Px = 0.05 and /?z = 0.35. For the former, we have 
used the topological ansatz with a string of operators at 
the virtual level, whereas the latter was obtained with an ex¬ 
citation that carries no string but transforms according to 
the odd representation of the virtual Z 2 symmetry (note that 
we did not impose this symmetry on the R-tensor in Eq. <§ 
explicitly). 


toric code, the flux of the excitations is encoded in the 
presence or absence of a string of operators, whereas 
the charge is encoded in the symmetry representation of 
the B tensoJ^. These topological characteristics are un¬ 
altered by the filtering. 

In Fig. we have plotted the elementary excitation 
spectrum in both the charge and flux sector for a spe¬ 
cific value of the filtering parameters. We can see that 
the filtering in one direction mostly affects one of the 
two excitations. This is reflected more clearly in Fig. 
where we have plotted the gap in both sectors along the 
Px = 0 axis. At the critical point = log(-\/2 -|- 1), one 
can observe that the gap to the one-flux and many-flux 
states closes, indicating the condensation of flux excita¬ 
tions. The charge sector remains gapped but will cease 
to exist beyond the transition, for reasons explained in 
Ref. HZ! 

Conclusions. This paper introduced the concept of 
post-PEPS methods as a systematic approach for study¬ 
ing the low-energy spectrum of two-dimensional quantum 
systems, as well as a new PEPS contraction scheme that 
enables to compute the necessary quantities. 

In the one-dimensional case, the local description of 
excitations in the MPS framework has led to the de¬ 
velopment of an effective particle picture of excitations 
on top of a strongly-correlated vacuuirPi]_ Qur methods 
open up the route to a two-dimensional generalization 
of this particle description of the low-energy dynamics 
and, more specifically, to the formulation of the scat¬ 
tering problem for effective particles in two dimensions. 
Applications such as multi-particle spectral properties, 
bound-state formation, magnon condensation and low- 
temperature behaviour seem within reach. 

We have also introduced a new contraction scheme, 
which will allow for PEPS algorithms based on the time- 
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FIG. 3. The lowest-lying excitation energies in the charge 
and flux sector as a function of /?z and with fix = 0, for 
which the model can be shown to have a phase transition 
at Pz = log(\/2 -I- 1) « O.sd^. The color represents trivial 
(red) or non-trivial (blue) flux of the excitations, whereas the 
symbol represents even (-I-) or odd (•) charge number. We can 
observe that the gap of the flux excitation (even charge in the 
non-trivial sector) closes at the phase transition, signalling the 
condensation of the flux excitations. We can even see the two- 
flux state (even charge in the trivial sector) coming down in 
energy close to the phase transition, although our variational 
ansatz is not suited for describing multi-particle states. The 
charge excitation (odd charge in the trivial sector) remains 
gapped. 


dependent variational principl^^ESl These will prove 
useful to simulate the real-time dynamics of quantum 
quenches, but might also result in improved algorithms 
for finding PEPS ground state approximations. Indeed, 
with our new effective environment we are able to com¬ 
pute both the gradient and the Hessian of the (global) 
energy functional, which can be used as an input for nu¬ 
merical optimization method^^. 
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SUPPLEMENTAL MATERIAL 


Transfer matrix and fixed points. In the first part of this Supplemental Material we study the fixed points of 
the linear transfer matrix T and corner transfer matrix C. For a translation-invariant MPS to satisfy the fixed point 
equation of the linear transfer matrix T 


r\j 

it is required that there exists a tensor (rectangle) such that 



holds up to small truncation errors. 

Let us now show that the “fixed point” of the corner transfer matrix C can indeed be approximated by introducing 
a single corner tensor S between the MPS tensors of the fixed points of respectively the horizontal and vertical linear 
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transfer matrix T; pictorially, we have 




Note, firstly, that the corner transfer matrix introduces two new sites at every application (similar to inhnite-size 
DMRG algorithm^2145ij therefore does not have a fixed point in the strict sense. Nevertheless, repeated application 
of C is likewise expected to result in a state with a converged (i.e. translation-invariant) structure, up to the corner 
itself. We therefore model the “fixed point” using the tensors of the fixed point of the linear transfer matrix and 
insert a new corner tensor. Applying the corner transfer matrix C once and using the tensor from Eq. gives rise to 



which shows that the original MPS tensors are indeed obtained after application of C, except on the two newly 
introduced sites. In principle, two new MPS tensors connected by a new corner matrix could appear. However, since 
the unique MPS tensor of the fixed point of T seem to capture the correct structure on the further sites, and these two 
can be assumed to have originated from previous applications of C, we can make the ansatz that these tensors should 
also be put on the two new sites. With this ansatz, we obtain a linear fixed point equation for the corner matrix 
itself, which corresponds to a simple eigenvalue equation. With the corner matrix as only variational parameters in 
the fixed point equation for C, we essentially have a linear subspace as ansatz and can therefore easily measure the 
error obtained by projecting onto this subspace. 

Excitation ansatz: some details. The variational ansatz for an elementary excitation is given by 


I (S)) =5] 

m,n 



The blue tensor H, containing all variational freedom in the ansatz, has the same dimensions as the ground state 
tensor A, yet the number of variational parameters is smaller because of a redundancy in the parametrization. Indeed, 
through simple insertion one can easily check that the tensors (where the green square is a, D x D matrix) 

Bo,x = □-HZI 

and 

Bo,y = ^ 

give rise to states with zero norm. Correspondingly, the effective norm matrix defined as 

(S')= 47r25(«, - n',)5{ny - <) B, 
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will have a number of zero eigenvalues (B is the vector containing all elements of the tensor B). In order for the 
variational subspace to be well-defined, we will always project out these null modes. In addition, since we want an 
excitation to be orthogonal to the ground state, we will also project out the component that is proportional to the 
ground state tensor A. 

Once we have defined the variational subspace, we can minimize the energy in order to find the best approximation 
to the true excitation. Since the subspace is linear, this can be done by solving the Rayleigh-Ritz problem 


H 


eff 


B = ujN 


Kj.Ky 

eff 


B 


where the effective Hamiltonian matrix 


is similarly defined as (Eq is the ground-state energy) 


(i?')| H-Eo- 4) B' B. 

The matrix elements of and contain two- and three-point functions that can be computed with a channel 

environment as explained in the main body. In the remainder of this Appendix we show how to do this in detail. 

Two-point functions. We will first look at the effective norm matrix, which yields the following double (infinite) 
sum 


n,n' 

Because the ground state is translation invariant, we can simplify to a single sum as 

= 47 r^i 5 ^(/?' — if) ^ QiKvm^-iKHn (m, 0) and B' at site (0,n')]. 

m,n' 

These terms are all two-point functions, which can be computed with our channel environment. One term looks like 



where the green (red) tensors indicate the locations of B and B'. In a first step we can contract the infinite channels 
by computing the fixed point of the “channel operators” and represent these as three-legged rectangles 
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The momentum superpositions can be worked out by inverting the relevant channel operator. For example, we can 
write one of the sums as 


giKv"! gj^g (m,0) and B' at site (0,0)] 


r 

n 

h 

u 

, r 


^2iKv *—*’ 






Jr 

-1 ^ gSiKv 


I 'I 



I 'i 




9-0-9 


I '1 


where 


i~r7 




-1 


(7) 


Applying this procedure to all geometric sums and keeping the local term amounts to 

I. ,1 


I . I 





■ r>VT^_n + 


I I 



I I 




-C 

D- 

-c 

)- 




-C 

3 — 




-c 


—c 

3- 


I . I 



I , I 



g-IKVg+iK/f 


-c 

3- 

-c 

3- 


q-IKv g-tliH 


r 


-r 

3 — 

4 



-c 

3- 

—c 

3- 
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( 8 ) 


where a green (red) square indicates the transfer operator with a B {B') on the ket (bra) level and the yellow square 
represents the situation where B and B' are on the same site; the big six-leg tensors represent inverses as in Q, the 
three-leg tensors are the fixed points of the channels and the diamond tensors are the corner matrices. 

The overlap of the Hamiltonian with respect to two excitations can be calculated easily for a frustration-free 
Hamiltonian: disconnected terms are automatically zero as the Hamiltonian annihilates the ground state locally. The 
overlap (<l>;;/[i3']| H reduces to a number of local contractions, which we will not write down (as they depend 

on the form of the Hamiltonian). 

Note that the non-frustration free case can be calculated equally well with our methods - provided we have a good 
PEPS representation of the ground state. Disconnected terms will generally be of the form 


I . I 



q 

p 




where the middle six-leg tensor represents a two-site Hamiltonian term squeezed between ground state PEPS tensors. 

Structure factor and sum rules in the AKLT model. Let us first study the static correlations of the AKLT 
state. In Refs. [H] and 05] it was argued that these correlations can be expressed as thermal averages of a related 
classical model in the same number of dimensions, contrary what one expects for generic quantum ground states. The 
real-space correlations then decay according to the 2D Ornstein-Zernike form 

s(n) = (4'(H)| |^(H)) oc e-l"l/« (n » 1) 

with 52 (t) the spin operator at site n and k* the oscillation period. In momentum space, the structure factor is thus 
given by 


3(^) 


s(K} = > e-*'‘'"s(n) 


= ('I'ol 5“;5o“(0) Ivko) 
1 


(9) 


for momenta close to i?* where s(k) reaches its maximum (we have defined the momentum spin operator 52 = 
Sn6*"^ "52). We can confirm this result with our methods. Firstly, we can calculate the correlation length by 
computing the gap of the linear transfer matrix T using the methods of Ref.^. We obtain the value ^aklt = 
— I/log|A^^| = 2.06491, in reasonable agreement with the value of « 0.52 in Ref. 071 In Fig. we compare 
the computation of the structure factor with our methods [Eq. (|^] with the form of Eq. © and observe very good 
agreement in a large portion of the Brillouin zone. 

The structure factor also shows up in the integrated spectral function. The latter is defined as 


5(^,0;)=^^ / dte' 




{^o\Sf{t)SS{0)\'So) 
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FIG. 4. The structure factor s{p) where p = — tt)^ + {ny — for the two-dimensional AKLT model along the axes 

Ky — TV (red) and = Kj, (blue) compared to the form in Eq. In the inset the (logjg of the) deviations are plotted. The 
plot shows that the classical Ornstein-Zernike form is accurate for a large portion of the Brillouin zone and that the structure 
factor is nicely isotropic around it* = (tt, tt). 


with S^{t) the spin operator at site n in the Heisenberg picture. By inserting a resolntion of the identity consisting 
of all excited states 1 = I 7 ) ( 7 I and only taking into account the one-particle states, we get the one-particle 

contribution to the spectral function 


5ip(«,..)= ^ |(7(«)|5o“|vk(H))|^ 

76 r(K) 

where r(K) is the set of all one-particle states | 7 (k)) with momentum k. Upon integrating the spectral function, we 
get the following sum rul^^ 

J = I ^ (vFol S%27r6{uj - H)S^{0) jd/o) 

= {'i>o\S%S^{0)\^o)=s{n). 






